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Abstract 

We use a recent 0(a 2 s ) calculation of the differential cross section for the 
production of heavy quarks in e + e~ annihilation to compute a few moments of 
the heavy quark single-inclusive production cross section. We verify that, con- 
trary to some recent claims, the leading and next-to-leading logarithmic terms 
in this cross section are correctly given by the standard NLO fragmentation 
function formalism for heavy quark production. 
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1 Introduction 



It is known that inclusive heavy quark production is a calculable process in pertur- 
bative QCD, since the heavy quark mass acts as a cut-off for the final state collinear 
singularities. Thus, the process 

e + e- ^Zh^Q + X , (1.1) 

where Q is the heavy quark and X is anything else, is calculable. Its cross section 
can be expressed as a power expansion in the strong coupling constant 

J oo 

-(*,£, m) = Y.a {n \x ) E ) m^)a n M , (1.2) 
ax n=0 

where E is the centre-of-mass energy, m is the mass of the heavy quark, p, is the 
renormalization scale, and 

As usual we define 

x = — — , (1.4) 

where q and p are the four-momenta of the intermediate virtual boson and of the 
final heavy quark Q. The cross section ( |1.2j ), normalized to the total cross section, is 
sometimes referred to as the heavy quark fragmentation function in e + e~ annihilation. 
When E/m is not too large, the truncation of eq. (|1.2|) at some fixed order in the 
coupling can be used to compute the cross section. On the other hand, if E ^> m, 
the n th order coefficient of the expansion will in general contain up to n powers of 
log (E/m), thereby spoiling the convergence of the expansion. These large logarithms 
can be resummed, according to the method described in ref. |1J. First of all, since 
log (E I m) is large, one is entitled to neglect the terms that are suppressed by powers 
of m/E. One then observes that, in this limit, the inclusive heavy quark cross section 
must satisfy the factorization theorem formula 

da ^ f 1 dz d&i -fx \ 

-(x, E,m) = Y: I — ^ (*> E, n) A * m) , (1.5) 



where d&i(z, E, p)/dz are the MS-subtracted partonic cross sections for producing 
the parton i, and D{(x, fi,m) are the MS fragmentation functions for the parton i 
into the heavy quark Q. In order for eq. ( pL.5|) to hold, it is essential that one uses 
a renormalization scheme where the heavy flavour is treated as a light one, like the 
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pure MS scheme. Thus d&i(z, E, \i)jdz has a perturbative expansion in terms of a s 
with nf flavours, where rif includes the heavy one. The scale \i is the factorization 
and renormalization scale. It should be chosen of the order of E, in order to avoid 
the appearance of large logarithms of E/fi in the partonic cross section. The MS 
fragmentation functions Z)j obey the Altarelli-Parisi evolution equations: 

— — (x,/i,m) = I — p ij (~>a a (ji)) Dj(z,n,m). (1.6) 
dlogfi 2 j Jx z \z J 

The Altarelli-Parisi splitting functions have the perturbative expansion 

P y (x,a a (/x)) = a a (jj)P^ (x) + a 2 (^P}}\x) + 0(a 3 s ) , (1.7) 

where p£° are given in ref. and P-p have been computed in refs. The only 

missing ingredients for the calculation of the inclusive cross section are the initial 
conditions for the MS fragmentation functions. These were obtained at the NLO level 
in ref. [|TJ by matching the 0(a s ) direct calculation of the process (i.e. formula ( |1.2| )) 
with the expansion of formula ( |1.5|) at order a s . They have the form 



b Q (x,no,m) = 5(1 — x) + a s (//o) dq\x, [Iq, m) + 0{a 2 s ) 
D g (x,ii ,m) = a s (no) d ( p(x,fi ,m) + 0(a 2 s ) , (1.8) 

all the other components being of order a 2 s . Thus, in order to compute the NLO 
resummed expansion, one takes the initial conditions eqs. ( |1.8| ), at a value of //o of 
order m, evolves them at the scale /j (taken to be of order E), and then applies 
formula Q1.5Q , using a NLO expression for the partonic cross section 

^{x,E,fi) = aj 0) (x) + al 1) (x,E,/i)a s (/i) + 0(a 2 ) . (1.9) 



For example, if the parton i is the heavy quark itself, one gets 

" n ~(x,E,fi) = 5(1 -x) + ag ) (x,£, y u)a s (/i) + 0(a 2 s ) , (1.10) 



fix- 
where we have normalized the cross section to 1 at zeroth order in the strong coupling 
constant. 

The procedure outlined above guarantees that all terms of the form (a s L) n (lead- 
ing order) and d s (a s L) n (next-to-leading order), where L is the large logarithm, are 
included correctly in the resummed formula. Observe that, at the NLO level, the 
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scale that appears in 6t s in eqs. ( |1.8|) and (|1.9|) could be changed by factors of order 
1, since this amounts to a correction of order a 2 . However, one cannot set fiQ = fi 
in eqs. (|1.8| ) (or // = m in formula Q1-9Q ), since this amounts to a correction of order 
a 2 L, and thus it would spoil the validity of the resummation formula at the NLO 
level. 

There is essentially no room for these considerations to fail. They are a conse- 
quence of the factorization theorem for fragmentation functions, which is quite well 
established [|, g. 

The validity of this procedure has however been questioned by Kniehl et al. in 
ref. ||. In their procedure, the heavy quark short distance cross section is replaced 
by 

^{x, E, M ) = 5(1 -x) + a s (E) [ag } {x, E, //) + d$(x, Mo , m)] + 0(a 2 ) , (1.11) 
and the initial condition by 

D' Q (x) = 5(1 - x) , (1.12) 

which is to be evolved from the scale no to the scale fi using the NLO MS evolution 
equations. This procedure differs at the NLO level from the standard procedure 
advocated in ref. PJ. The difference starts to show up in the terms of order a 2 L. 
Recently, we have completed a calculation of heavy quark inclusive production at 
order a 2 [10]. Using this calculation, we are in a position to verify explicitly the 



approach of ref. and thereby dismiss the approach of ref. In the next section 
we describe the procedure we followed in detail. 



2 Calculation 

Instead of dealing with the realistic case of Z/7 decay, we perform the calculation 
for a hypothetical vector boson V that couples only to the heavy quark with vectorial 
coupling. 

We introduce the following notation for the Mellin transform: 

f(N)= f 1 dxx N - 1 f(x). (2.1) 
Jo 

We adopt the convention that, when iV appears instead of x as the argument of 
a function, we are actually referring to the Mellin transform of the function. This 
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notation is somewhat improper, but it should not generate confusion in the following, 
since we will be working only with Mellin transforms. The Mellin transform of the 
factorization formula (|1.5| ) is given by 



a(N, E } m) = J2 °i(N, E, y) D^N, fi, to) , (2.2) 

i 

where i 

a(N,E,m) = f 1 dxx N - 1 ^-(x,E,m) , (2.3) 
Jo dx 

and a similar one for &%(N, E, fi), and the Mellin transform of the Altarelli-Parisi 
evolution equation (|1.6|) is 



^dhg^ = 5>'M [ P « 0) W + + 0(a 2 s )] bj(N,fx,m) . (2.4) 

We want to obtain an expression for a(N, E, m) valid at the second order in a s . Thus, 
we need the solution of eq. Q2.4|), with initial condition at // = yo, accurate at order 
a 2 . This is easily done by rewriting eq. fl2.4|) as an integral equation 



Di(N, /i, m) = Di(N, /i , to) 

+ £ T^log/i' 2 a s (/V) \P^(N) + P^(N)a s (fi')} D^N^'.m) . (2.5) 

• J Un 



The terms proportional to a 2 can be evaluated at any scale (/i or /i ), the difference 

^,3 



being of order a 6 s . Factors involving a single power of a s can instead be expressed in 
terms of a s (/xo) using the renormalization group equation 

Mo 



a 8 (ji') = a a (jJLo) - 2vr b a 2 (/j ) log -y + 0(c^(/i O )) 



6o = "0,-4,^ (2.6) 

1Z7T 

where n/ is the number of flavours including the heavy one. Equation ( |2.5| ) then 
becomes 

Di(N,ti,m) = D i (N, f i ,m)+J2 I* dlogy' 2 a s {y Q )P§\N)b j {N,y\m) 

+ E« s 2 W4¥)A{^ft,m) log 4 
j Mo 

- 27r6oE«sM4 0) (^)A'(^Mo,m)^log 2 4- (2-7) 

z Mo 



-6- 



We now need to express Dj(N, /i' ,m) on the right-hand side of the above equation 
as a function of the initial condition, with an accuracy of order a s . This is simply 
done by iterating the above equation once, keeping only the first two terms on the 
right-hand side. Our final result is then 

u 2 

Di(N, fi, m) = Di(N, /i 07 m)+J2 aa(jio)P§\N)Dj(N, /i , m) log ^ 
+ E <?Mp£\n) P$\N)D 3 iN, aio, m) \ log 2 ^ 

kj 1 

+ J2^o)PiP(N)D j (N^ ,m) log ^ 

- 2nb Y,^o)P{j\N)D J (N^ 0} m)l\og 2 ^. (2.8) 

Since the initial condition is 

D 3 (N,fi ,m) = 5 jQ + a s (fi )df\N,fi ,m) + O(a 2 s (fi )), (2.9) 
eq. (|2.8|) becomes, with the required accuracy: 

i o z^ 2 

Di(N,n,m) = 5 iQ + a a (jio) df> (N, fi ,m) + a 3 (fi )P$(N) log — 

A*o 

+ £^(,io)pW(JV)^(JV )ttl ,m) log ^ 

+ E« s 2 M4 0) (iv)p(g(iv)iiog 2 4 

+ ^(/io^AOlog^ 

- 27r6 « 2 MP J S ) (iV)ilog 2 4- (2-10) 

Re-expressing a s (/io) in terms of /x, we get 

i i 

+ a.(/*)P$(JV) hg^ 2 +J2^)P?(N)df\N,^m) log ^ 

^0 1 
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/'o 



/'o 



f2.11^ 



The partonic cross sections are given by 

&i(N, E, /i) = 5 lQ + % + a a (ji) d^(N, E, //) + C(a s 2 (//)) , (2.12) 

where a!p vanishes unless i is either Q, Q or g. Thus, combining eq. fl2.12|) with 
eq. ( |2.11| ) according to eq. ( |2.2j ), we obtain 



a{N,E,m) 



1 + a s (n) 



a$(N, E, y) + d$(N, y ,m) + P™(N) log , 



7(1)/ 



>(0), 



+ a 2 (/i)<j J2#\N,E,n)P$ log^ + 2Trb Q d$(N,n ,m) log^ 



,(0) 



7(1)/ 



/'0 



/'o 



+ EK'W+^W] 4\ N >*>,m) log 



+ 



fcQ' 

,2 



/'0 



The above formula should accurately describe the terms of order a s L, a s , a 2 L 2 and 
a 2 L. Terms of order a 2 , without logarithmic enhancement, are not accurately given 
by the fragmentation formalism at NLO level, and have consistently been neglected. 

The lowest order splitting functions are given by 



(2.13) 



P^(N) = C F 



2S 1 (N) 



P®(N) = C F 



2 iv(jv + i; 

2 + N + N 2 



P°MN) = Tp 



N(N 2 - 1] 
2 + N + N 2 



N(N + 1)(N + 2) 
where, restricting ourselves to integer values of N, 



N 



j=l J 



(2.14) 



(2.15) 
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pW(iV) and ?W (iV) are given by 



^S(iV) = + P g > g (N) , PghN) = PfJN) + P q , g (N) , (2.16) 



where the non-singlet components are given by 

P$%(N) = Pqq(N) + P QQ( N ) + Pqq( N ) > P f Q ( N ) = P ^ N ) + P^ N ) > 

= [Pf(N) + A(N)} , P^(iV) = \C F C A P G {N) , 

Pq Q (N) = n f C F T F P NF (N) , 



pg(iV) = C F P A (N) 



p^( n ) = -\c f c a p a {n) 



(2.17) 



and 



D , , n „ „ 8 + 44iV + 46iV 2 + 21iV 3 + 14iV 4 + 15iV 5 + 1(W 6 + 2N 7 
Pq'q{N) — —C F T F N3(N + mN + 2y(N-l) ' (2 ' 18) 

P F (N), A(N), Pq{N) and Pnf(N) were taken from the appendix of ref. U and 
P A (N) is given in eq. (5.39) of ref. [Q. We have obtained our explicit expression for 
Pq'q(N) using the equation 



P- 



Pq'q 



pNS 
QQ 



pNS 
QQ 



2ra 



(2.19) 



where Pqq is the singlet component] 2 ], calculated in ref. ||. Equation (|2.19[) is easily 
seen to follow from eqs. Q2.16 ) and from eqs. (2.42) of ref. |J. 

The expressions for Oq* and Sq are respectively given in eq. (A. 12) and (A. 13) 
of ref. []]]. The coefficient a^> can be obtained by performing the Mellin transform of 
the expression cx, g + CL, g , where or, g and CL, g are given in eq. (2.16) of ref. f9j. Thus 



-2 



V 2 



log — + 4 



' 2(2 + iV + iV 2 
N(N 2 - 1) 

^ lSl (7V-l)-| Sl (7V) + — 



+ 



(iV-1) 2 N 2 (iV + 1) 2 
1 



Si(N+l) 



d q 1] (N, 



P^(N) log 



m 2 



[2.20) 



2 We warn the reader that, sometimes, in the literature, the notation P s is used for the "sea" 
component, and Pqq is used for the singlet one. Here we use Pqq for the full QQ splitting function. 
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In order to make a more detailed comparison with our fixed order calculation, we 
separate the 0(d 2 ) contributions to a(N,E,m) according to their colour factors. 
Choosing for simplicity \i = E and /iq = m, and using the notation 



we write 



d Q/g( N ) = d Q ] /g( N > m )L=m , 



(2.21) 



1 + a 8 (E) A(N, E, m) + aj(E) B(N } E, m) 

E 2 



7(1)/ 



a(N,E,m) 

A(N,E,m) = clq\n) + dQ } (N) + Pqq(N) lo. 
B(N, E, m) 

B CF (N,E,m) 



5 CF (iV, P, m) + Pc A (iV, m) + B n . (N, E, m) + B Tp (N, E, m) 



[P { Q ] Q (N) [d$(N) + a$(N)} + P&(N) + Pfj(iV)} log 



P 2 



1 r 



p 2 



B CA (N,E,m) 



P&(N) + Pg(iV) + ^C A d$(N) 



11 
~6~ 



log^ 



12 



B n .(N,E,m) 



P% Q {N)-\n f T F d%{N) 



log 



P 2 



B TF (N,E,m) 



i(i)/ 



(AO Pig + 2 P$(JV) 4 X) (iV) + 2 P, q {N) log 



,(0) 



P 2 



(0), 



P 2 



f2.22) 



where the subscripts CV, CU, rif and Tp denote the C F , C f Ca, rifC F T F and C F T F 
colour components. 

The fixed order calculation of ref. JTO] can be used to compute the cross section 
for the production of a heavy quark pair plus one or two more partons, at order a 2 . 
We separate contributions in which four heavy quarks are present in the final state, 
from those where a single QQ pair is present together with one or two light partons. 
These last contributions were computed only in a three-jet configuration, and they 
are singular in the two-jet limit, that is to say, when x—*l. Furthermore, the virtual 
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corrections to the two-body process V — > Q + Q are not included in our calculation. 
In order to remedy for these problems, we proceed as follows. The 0(a 2 ) inclusive 
cross section for V — > Q + Q + X, can be written symbolically in the following form 

^- = a ^5(l-x) + a s j dY a {l \x,Y) + a 2 s J dY a\ 2 \x,Y) + 2 J dYa£\x,Y) 

(2.23) 

where Y denotes all the other kinematical variables, besides x, upon which the final 
state may depend. We assume fj, — E, and we do not indicate, for ease of notation, 

(2) 

the dependence upon E and m of the various quantities. The term a\ arises from 

— (2) 

final states with a single QQ pair plus at most two light partons, while a h arises 

— (2) 

from final states with two QQ pairs. The factor of 2 in front of the a h contribution 
accounts for the fact that we may detect either one of the two heavy quarks. The 
moments of the inclusive cross section can be written in the following way 



<j(N) 



jv_i da 
dx 



,/,,,- - J -_ = a + a s JdxdY(x N ~ l - l)a (1 \x,Y) 



J dxdY(x N ~ l -l)a?\x,Y) +2 J dx dY {x N ^ - \) af^x.Y) 



(2.24) 



where 



a = a (0) + a s Jdx dY a (1) (x, Y) + a 2 s J dx dY a{ 2) (x, Y) + J dx dY a% ] (x, Y) 

(2.25) 

The expression for cr(N) can now be easily computed with our program, since the 
(x N ~ 1 — 1) factors regularize the singularities in the two-jet limit, and suppress the 
two-body V — > Q + Q virtual terms. Furthermore, in the massless limit 



o = \. + 2rx„{E) + c<X 2 H {E) 0[ T -^\-0(ai) 



2\ 



(2.26) 



where c is a constant. In fact, the 0(a 2 s ) term does not contain any large logarithm, as 
long as a s is the coupling with rif flavours, including the heavy one 3 . Reintroducing 
the energy and mass dependence, we have 



a(N,E,m) 



a s (E) L(N, E, m) + a 2 s (E) M(N, E, m) 
+ ca 2 (E) + o(^)+0(al) 



3 If instead the cross section formulae are expressed in terms of a{ nf ^ 

-,( n f — 1) z' tp\ , 4rp_,„„ B 2 -2 , frs(-0. 



we have 



a = l + 2os 



>(E) + lT F l0g^aj + O(aj). 
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L(N, E, m) = 2 + j dxdY (x*- 1 -l)a {1) (x,Y) 
M(N, E,m) = J dx dY - l) af ] (x, Y) 

+ 2 J dx dY (a;^ 1 - i) off (x, Y) . (2.27) 

We have calculated L(N, E,m) and M(N,E,m) numerically, using E = 100 GeV 
and m — 8, 4, 3, 2.5, 2, 1.5, 1, 0.6, 0.5, 0.4, 0.2 GeV, for a vector current coupled 
to the heavy quark. We expect that, for small masses, A(N,E,m) should coincide 
with L(N, E,m), and M(N, E,m) should differ from B(N,E,m) by a mass and 
energy independent quantity, since such term is actually beyond the next-to-leading 
logarithmic approximation. We find very good agreement between A(N, E, m) and 
L(N, E, m). We present the results for M(N, E, m) separated into the different colour 
components 

M(N, E, m) = M Cf (N, E, m) + M Ca (N, E, m) + M nj (N, E, m) + M Tp (N, E, m) . 

(2.28) 

In Fig. [I] we have plotted our results for Mq a (crosses with error bars) and for Bq a 
(solid lines). An arbitrary (iV-dependent) constant has been added to the curves for 
o 

-1000 

■+-> 

PI 

0) 

o -2000 

o 
o 
o 

ia -3000 

-4000 



4 6 8 10 12 

log E /m 

Figure 1: CfCa component of the a 2 s coefficient in a(N, E,m), as a function 
of log E 2 /m 2 , for N =2, 5, 8 and 11. 

Bq a , in order to make them coincide with the numerical result for m/E = 0.015. We 
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find, as the mass gets smaller, satisfactory agreement for all moments. Notice that, 
as intuitive reasoning would suggest, for higher moments we need smaller masses to 
approach the massless limit. In Figs. |2|, [3] and f| we report the analogous results for 
the remaining colour combinations. Again, we find satisfactory agreement. 
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Figure 2: Same as in Fig. for the nfCpTp component. 



If one follows the procedure proposed by Kniehl et al., eqs. (|2.22|) are modified in 
the CfCa and in the rifCpTp coefficients. More specifically, the terms proportional 
to dg^ all disappear from the expressions of the CfCa and of the npCpTp coefficients. 
In fact , by inspecting formulae ( 1.11 ) and ( |1 . 1 2j ) , and the derivation of eq. (|2.13|) , we 
see that the only relevant difference between the two approaches is that the term 
a s (//o) do is replaced by a s (/i) dq\ which, using the renormalization group equation, 
amounts to a difference of — 2tt b o? s dg log \i 2 j precisely what is needed to cancel 
the term of the same form appearing in eq. (|2.13| ). The modified result of Kniehl et 
al. is also shown in Figs. |I] and |2| (dashed lines). It is quite clear that their approach 
does not work. 



13 



2500 



2000 



Solid: standard formalism 
X: numerical result 



o 1500 



o 
o 
o 

Nm 1000 

its 



500 




Figure 3: Same as in Fig. [J, for the C F component. 
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Figure 4: Same as in Fig. for the CfTf component. 



3 Conclusions 



In the present work, we have verified at order a 2 s the NLO fragmentation function 
approach to the computation of the heavy quark fragmentation function given in 
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ref. 0. Besides excluding an alternative that has been proposed in the literature §|, 
we have given a verification of several ingredients that go into the fragmentation 
function formalism, such as the initial conditions, computed in ref. [l|], the NLO 
splitting functions in the time-like region H, and finally we have also performed a 



further test of the validity of the fixed order calculation of ref. fllO 
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